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Abstract 

We study theoretically the sound propagation in charge- and spin-density 
waves in the hydrodynamic regime. First, making use of the method of co- 
moving frame, we construct the stress tensor appropriate for quasi-one di- 
mensional systems within tight-binding approximation. Taking into account 
the screening effect of the long-range Coulomb interaction, we find that the 
increase of the sound velocity below the critical temperature is about two 
orders of magnitude less for longitudinal sound than for transverse one. It is 
shown that only the transverse sound wave with displacement vector parallel 
to the chain direction couples to the phason of the density wave, therefore we 

expect significant electromechanical effect only in this case. 
72.15.Nj, 71.45.Lr, 75.30.Fv 
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I. INTRODUCTION 



Several aspects of the collective transport associated with the phason in charge- or spin- 
density waves (CDW or SDW) are still not well understood. One of the intriguing phenomena 
is the electromechanical effect observed in CDW |l]-[5| and more recently in SDW ||. First 
of all, most of the elastic moduli increase upon entrance into the CDW or SDW state, often 
with a sharp dip at T c , the transition temperature. Second, some of the elastic moduli in 
CDW or SDW soften when the density wave is depinned by an external electric field in 
excess of the depinning threshold field E?- Third, the change in the elastic moduli due to 
depinning of the density wave depends on the frequency u> of the flexural vibration f7j, and 
decreases like uo~ p with This behavior is similar to the frequency dependence of the 

change in the dielectric constant upon depinning in CDW and SDW ||. 

We have shown earlier || in the collisionless limit that the hardening of the elastic 
constants can be understood in terms of the reduction in the quasiparticle screening of 
the ion potential due to the formation of the density wave state. The electromechanical 
effect was interpreted as an additional screening contribution from the collective mode of 
the density wave condensate (phason) liberated by depinning. A later extension of that 
theory to the experimentally more relevant hydrodynamic limit [[KJ did not modify the 
above picture qualitatively. However in these papers it was assumed that the phonon simply 
couples to the electronic density, and the effect of the long-range Coulomb interaction was 
neglected. 

The purpose of the present paper is twofold. First, we shall develop the theory of 
the electron-phonon coupling for a strongly anisotropic system, which will enable us to 
distinguish between the behavior of transverse and longitudinal sound waves propagating in 
various directions. Second, we shall include the effect of the long-range Coulomb interaction 



following Kadanoff and Falko [p| . In Section II we construct the electronic stress tensor 
(which couples to the deformation tensor of the sound wave) for a quasi-one dimensional 



system following the method of comoving frame In Section III we concentrate on 
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the quasiparticle contribution to the stress tensor correlation functions corresponding to 
the pinned case. Section IV is devoted to the examination of the coupling of the stress 
tensor to the phason, which is relevant to the electromechanical effect. Our conclusions are 
summarized in Section V. A preliminary report on this work has already been published 
elsewhere ]TB[]. 



II. ELECTRON-PHONON COUPLING 



Following Tsuneto [12| let us assume that both the ionic potential and the electronic 
wavefunction are deformed in the slowly varying sound field u(r, t) = u cos(qr — ujt) imposed 
externally (extreme tight-binding limit): 

K(r)-> V[r-u(r,t)\ 

tjj(r) ->^[r-u(r,t)](l + Vu)" 1 / 2 . (1) 

This displacement is generated by the unitary operator U = exp[— u(r, £)V]. Therefore the 
deformed wavefunction is an eigenfunction of the transformed Hamiltonian h = Uh U + , 
where ho = e(— zV) with e(p) being the zone-periodic electronic energy spectrum. The 
transformed Hamiltonian h is then expanded in terms of the deformation tensor Vjttj, which 
is much smaller than one, even though the displacement u itself may be many times of the 
lattice constant for sound propagation. We obtain h = ho + h e i- p h, where the Hamiltonian 
for the electron-phonon coupling is given by 

hei-ph = ^2(^iUj)iVjVi(-iV), (2) 

with v(p) = de(p)/dp, the velocity of the Bloch electron. The matrix element of h e i_ p h 
between Bloch states is evaluated as 

(p + q\h el _ ph \p) = -iJ^QjUiTijip), (3) 

where the stress tensor 
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Tij(p) = mvi(p)vj(p), (4) 
and m is the bare electron mass. This expression generalizes the stress tensor used for an 



isotropic metal [11 



In orthorombic symmetry the sound wave polarized in the i direction and propagating 
in the j direction couples to the Ty component of the stress tensor, and in order to deter- 
mine the effect of that coupling on the frequency (or velocity) of the sound, we have to 
evaluate the appropriate stress tensor correlation function ([ry, t^]}. Once this is known, 
the renormalized sound velocity can be calculated in the weak coupling limit as 

c = c {l-([r,r])/2Mc 2 }, (5) 

where Cq is the sound velocity without electron-phonon coupling, M is the ion mass, and 
for clarity we have suppressed the indeces both for the stress tensor component and for the 
sound velocity 

For a highly anisotropic (t a t& t c ) tight-binding dispersion 

e(p) = -2t a cos(ap x ) - 2t b cos(bp y ) - 2t c cos(cp z ) - //, (6) 



widely used for CDW and SDW materials fl4| , the velocity and stress tensor components are 



easily obtained, and their values on the open Fermi surface can conveniently be expressed by 
the component of p perpendicular to the chains (x direction). However, since the Green's 
functions for CDW and SDW are usually written in the left-right spinor representation 
flQ| , involving measuring momenta from ±Q/2 with the density wave wavevector Q = 
(2pF,7[/b,7[/c), we should express the stress tensor elements in a compatible manner. It 
turns out, that for each stress tensor component the term proportional to the unit matrix 
dominates, therefore we have 

t xx = mv 2 F {l - [t b / 't a sin(ap F )} 2 sin 2 {bp y )} 
r yy = mVy[l + cos(2bp y )] 

r xy = mv F v y y/2cos(bpy). (7) 



Here pf is the Fermi momentum, vf = 2at a sin(apF) is the Fermi velocity in the chain 
direction, while v y = \f F ibt\ ) <C Vp is a typical velocity in the perpendicular direction. In the 
followings we restrict our study to the x — y plane, since behavior involving the z direction 
should be similar to that of the y direction. We note, that all components of the stress 
tensor depend on momentum through the combination if = bp y only. 

III. PINNED DENSITY WAVES 

In this section we consider sound propagation with no applied electric field. The density 
wave is pinned, therefore the condensate is unable to contribute to correlation functions, 
including that for the stress tensor. Mathematically this situation can be simulated by 
setting the coupling of the stress tensor (and of the density) to the phason to zero. Then the 
stress tensor couples only to the density fluctuations, resulting in the well known Coulomb 
screening [|TTj| : 

, w j» iw -A\ ([r, n])o([n,r]) 

(t,t) = (t,t) yj- — YTT\ n~- ^ 

q z /Ane 2 + {[n, n\) 

Here n stands for the electronic particle density, and ([A, B]) denotes a correlation function, 
in which only the effect of impurity scattering is taken into account. 

The density correlator ([n,n]) in the presence of impurity scattering was evaluated in 
|T0| . A straightforward extension of that calculation confirms that under the circumstances 
of the sound experiment (Iq <ti 1, where / is the mean free path) the stress tensor correlator 
has two distinct contributions: 

(Mo = (T(<p))l([n,n]) + ([5r( V )}\([n,n])r ertex . (9) 

The first contribution features only the average of the stress tensor (t((^)) v , = 
(27r) _1 JZ n dipr (if), and is proportional to the diffusive (vertex corrected) density correla- 
tor. The second contribution containes only the fluctuating part of the stress tensor com- 
ponent 6r(<p) = r(ip) — (t((p))<p, and therefore it is proportional to the density correlator 
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([n,n])Q° vertex calculated without vertex corrections. According to the same argument, only 
the average of the stress tensor couples to the density, therefore 

([r,n}} = ([n,r}} = (r(if)) v ([n,n}} . (10) 

Combining Eqs.(8)-(10) we see that in the long wavelength limit appropriate for the sound 
experiment (q ~ 1/L, where L is the sample size), the average part of the stress tensor 
(s-wave component, proportional to density) is completely screened out by the Coulomb 
interaction, and only the fluctuating part contributes to the correlation function: 

([r,T}) = ([5r( V )] 2 ) v ([n,n})r ertex . (11) 

The situation here is the same as in the electronic Raman scattering, where the long-range 
Coulomb interaction suppresses the density (charge) fluctuations, and only non s-wave chan- 



nels survive |[15|| . 



The evaluation of ([n, n ^™ overtex can be done starting with the results of [10|. The 
calculation is rather technical, therefore we delegate it to the Appendix, and we give here the 
results only. Without vertex corrections there is no diffusion pole, and both the wavenumber 
q and the frequency ui could be set to zero. However we keep a finite (but small) frequency 
for finite imaginary part of the correlator: 



([n,n])™ vertex = N F lTqpil - f \ (12) 

LO + lT qp 

where Np is the density of states at the Fermi surface. The corresponding "unrenormalized" 
condensate density / (for the general formula see Eq.(32) in the Appendix) and quasiparticle 
damping T qp are evaluated in two limiting cases, close to T c and close to zero temperature 

as 

f(T - T c ) = - 2(^)V'(i + 5&0 « 7 C(3)(^f ) 2 

f(T -»• 0) = 1 - 37ra/i6, (13) 

and 



f qp (T -> T c ) = 2r 

f fT (V) = 9lT Aq8/3 G r G/T fUl 

1 gpl J — »■ — 3 2 u 2( 5 _ 4a 2/3) T e • l 1 ^ 

Here A is the density wave order parameter, r = rV + rs/2isa combination of the impurity 
forward and backscattering rate, a = T/A, = 1 — a 2 / 3 and G = Auq is the density wave 
gap. 

As it is seen from Eq.(13), / increases linearly in (T c — T) below T c , but it is slightly less 
than one at T = (T is usually an order of magnitude smaller than T c ). Nevertheless, the 
temperature dependence of the sound velocity in the pinned case will still be qualitatively 
the same as in the collisionless limit || . The relative change of the sound velocity compared 
to the normal state (c n ) is easily obtained from Eqs.(5), (11) and (12) as 

(c - c n )/c = Xf, (15) 

where the effective coupling 

A = ik {[Srm% - (16) 

Using Eq.(7), these effective couplings for the various sound waves are 

^ _ N F {mv 2 F ) 2 r tb 14 



A 



16AfCp L t a sin(app) ' 

Np(rriVy)' 2 
VV 4Mc% 



\ _ N F (mv F v y ) 2 / 17 \ 

A *y ~ 2M? • \ L{ ) 



Since v v /vf ~ h/t a m 1/10 in many quasi one dimensional materials, we expect that the 
relative increase of the sound velocity below T c will be a factor 10 2 smaller for longitudinal 
sound than for transverse sound. 



IV. ELECTROMECHANICAL EFFECT 



If an external electric field in excess of the threshold field Et of the nonlinear conductivity 
is applied in the chain direction, then the condensate is depinned and is able to contribute 
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to various correlation functions ||. The best known example is of course the conductivity 
itself, but the situation is the same for the stress tensor correlator as well. The collective 
contribution to ([r, r]) can be obtained if we allow both the stress tensor and the density to 
couple to the phason (in the previous section this coupling was blocked due to pinning). In 
this case the stress tensor correlator has another contribution ([t,t]) co11 in addition to the 
one calculated in the previous section, namely: 



xColl _ U( [r,5A]) c ^({5A,r ]) c ^ 
1 - U{[5A,5A]Y 



I r_ _i \ coil 

(l T ' r ]> = 1 rr/r.A rAlVCoui • ( 18 ) 



Here U is the on-site Coulomb repulsion responsible for the formation of the SDW state 
(for CDW it should be replaced by the phonon propagator, but that does not affect our 
conclusions), 6 A is the phase fluctuation of the order parameter, and ([A,B]) is the 
correlation function of quantities A and B including the effect of the long-range Coulomb 
interaction (like in Eq.(8)). As we have seen earlier, in the long wavelength limit this yields: 

{[AB]) Coui = {[AB])o _ ([A,n])o{[n,B}) /{[n,n})o. (19) 

First we consider if the allowed coupling to the phason actually takes place for various 
sound waves. According to Eq.(18) we need to examine ([r, 5A]) Coul , which is given by 
Eqs.(19) and (10) as 

([r,6A}) Coul = (MA]) - <t(<^(MA]> . (20) 

The density-phason correlator ([n, 5A])o was evaluated in [|K|]. Here we only reiterate that 
result in the limit of experimental interest Iq -C cq/vf (dynamic limit): 

([n,6A]) = i N F { -^^f d , (21) 

where in our two dimensional geometry the wavenumber q appears in ((p) = vfQ x + 
q y cos(f, and the condensate density in the dynamic limit is given by: 

f,(T -> T ) = a2 ib'(- + —) w 7rA " 

fd(T - 0) = 1. (22) 



We recall that the same fd appears in the current-phason correlator ([j,5A}) , as well as 



in the phason propagator ||16|| . Note that fd increases from zero much faster below T c than 
/ does, and that at zero temperature it saturates exactly to 1. The stress tensor-phason 
correlator can be calculated similarly: 

(|r, dA\) = iN F — f d . (23) 

Now we shall examine Eq.(20) for different sound waves in order to determine if there is 
a collective contribution to the corresponding stress tensor correlator. We consider longitu- 
dinal and transverse sound waves propagating in the x and y directions. Clearly, the second 
(screening) term in Eq.(20) is nonzero only for the longitudinal sound propagating in the 
chain direction (q || u || x), in which case it completely cancels the first term, leading to 
no collective contribution. The other longitudinal sound propagating perpendicular to the 
chains (q || u || y) does not couple to the phason either, because (T yy (ip) cosy?)^ = (see 
Eq.(7)). The coupling of the transverse wave propagating in the chain direction (q || x and 
u || y) is controlled by (r^^ip))^ = 0, yielding again no collective contribution. This means 
that in all of the above three cases there will be no electromechanical effect. 

For the rest of this section we will concentrate on the only interesting case, when the 
transverse sound propagates perpendicular to the chains (q || y and u || x). In this case 
there will be coupling to the phason, since 

([r xy , 5A}) Coul = ([r xy , 5A}) = t N F ^mv F v 2 y q y (24) 

is nonzero. Now we have to consider the denominator in Eq.(18). Since ([n, SA])q = for 
q x = 0, therefore ([5 A, 5A]) Coul = ([SA, SA]) , and we can use the result for the phason 
propagator calculated in pi, which in our case reduces to 



1 - U([5A,5A]) = ^0\M 2 ~ iu>T ph ). (25) 
Here T p h is the phason damping rate and is given by 

T ph (T -> T c ) = 2Y B 

W ^ °) = SKf e - G/T (26) 
9 



The phason damping freezes out for low temperature, and approaches 2Tb at T c . Note 
the discontinuity in Tp^ at T c (approaching from above Y p h « 27r 3 T/7£(3)), which is the 
consequence of the finite order parameter A below T c exceeding almost immediately the 
energy scale set by u and vq. 

We are now able to write down the total correlation function for this sound wave in the 
unpinned case. Using Eqs.(ll),(18), (24) and (25) we obtain 



The above equation means that the collective contribution (the second term on the right 
hand side) due to the moving condensate recovers some of the screening of the ion motion 
lost because of the decrease in the number of quasi particles. In fact, at zero temperature it 
overcompensates somewhat, since fa > f. Therefore we expect the electromechanical effect 
for the transverse sound polarized in the chain direction only. According to Eq.(27), close 
to T c the electromechanical effect on the sound velocity should be much smaller than the 
temperature effect, while at low temperatures the softening is somewhat bigger than the 
hardening was upon cooling. Although the collective contribution in Eq.(27) does have a 
frequency dependence, it does not appear to describe the suppresion of the electromechanical 
effect when the frequency is increased [0 . This is rather puzzling, although the above results 
for sound propagation may not translate literally to the flexural experiment. 



We have derived for the first time the appropriate stress tensor for quasi-one dimen- 
sional electron systems. Under conditions of a sound velocity measurement (Iq <C 1) the 
long-range Coulomb interaction has a simple role to suppress the s-wave channel as in the 
theory of electronic Raman scattering. In highly anisotropic systems like Bechgaard salts 
the increase of the sound velocity in the density wave state is two orders of magnitude 




x 



N F (mvFVy) 2 x 

»fq P (l-/) _i_ (■"yQyffd 



(27) 



w+if w (v y q y ) 2 -iuT ph 



V. CONCLUSIONS 



10 



smaller for longitudinal sound waves than for transverse ones. We also find, that a sound 
wave with polarization perpendicular to the chain direction can not couple to the phason, 
because the density wave condensate can only move parallel to the chains. The coupling of 
the longitudinal sound propagating in the chain direction is screened away by the Coulomb 
interaction, which leaves us only the transverse sound wave propagating perpendicular to 
the chains as the one which does couple to the phason, and shows the electromechanical 
effect. 

Recently Britel et.al. measured the elastic constant C44 of (Tai-iNb^Se^I at 15MHz in 
the geometry u || x |17|], and found a relative reduction of order 10 -4 in the presence of an 



electric field approximately IOEt- This seems to be consistent with our analysis, although 
the observed effect appears to be a little too small. 
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APPENDIX 

We evaluate here the density correlation function ([n, n ^™ overtex (Eqs.(12)-(14) in the 



text). We start with the corresponding thermal product (See [|iOf ): 



([n,n})r ertex = N F [1 - 7rT^F(iu n ,iu' n )}, (28) 

n 

where Np is the density of states, and u n and u' n are related to the Matsubara frequencies 
u n and J n = io n ^ v by 
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uj n /A = u n [l-a(u 2 n + l)- 1 / 2 } 



(29) 



with A the order parameter, a = T/A and T = Tp + r#/2 a combination of the impurity 
forward and backscattering rates. Neglecting vertex corrections in the relevant formulas in 
TOfl leads to: 



F(u,u') 



i I l+uu' 

L (l-U 2 ) 1 /2(l- U '2)l/2 



(30) 



A[(l-U 2 ) 1 /2 + (l_ n '2)l/2]' 

where u and u' are analytic continuations of iu n and iu' n , and in the absence of the diffusion 
pole the wavenumber q was already set to zero. 



While evaluating the correlation function we follow the standard method [18]. Expanding 
up to linear order in u we obtain 



([n, n})^^ = N F (1 - f + iul), 



(31) 



where 



and 



/ = ^E(^ + ir 3/2 , 



(32) 



x 



x 



j _ J_ roodE ,-2 (E_ \ 
1 ~ 2A JG 2T CObI1 \2TJ 

- Refl - u 2 )- 3 / 2 



(33) 



Here ti = (1/2)[1 + (|m| 2 + 1)/|m 2 - 1|], and G = Au 3 is the gap with u 2 = 1 - a 2 ' 3 . The 
above equations can be evaluated in two limiting cases, with temperature close to T c and 
close to zero, and can be brought to the form of Eqs.(12)-(14) in the text. 
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